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Abstract 

We show how the longstanding problem of the collapse of the charge-exchange QRPA 
near the physical value of the force strength can be circumvented. This is done by in- 
cluding the effect of ground state correlations into the QRPA equations of motion. The 
corresponding formalism, called renormalized QRPA, is briefly outlined and its conse- 
quences are discussed in the framework of a schematic model for the two-neutrino double 
beta decay in the ^^^Mo ^^^Ru system. The question of the conservation of the Ikeda 
sum rule is also addressed within the new formalism. 
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Double beta {(3f3) decays occur in medium-mass nuclei that are rather far from closed 
shells. The nuclear structure method most widely used in the evaluation oi (3(3 rates for 
two-neutrino decay mode {(3(32v) as well as for the neutrinoless mode {(3(3qu) is therefore the 
quasiparticle random phase approximation (QRPA) These calculations, in which the 
(3(32u matrix elements J^2y are approximated by their J'^ = l"*" component {i.e., J^2u — 
M.2u{J^ = 1^))) explain the smallness of the measured transition rates. [| However, 
the actual value of J^2v depends sensitively on the strength g^^ of the particle-particle 
force in the 5* = 1, T = channel. For realistic forces of finite range M.2v passes 
through zero near g^^ = 1 i.e., near the physical value of this coupling constant. This 
feature makes the actual value of Ai2iy rather uncertain. What is still more distressing, 
is that QRPA collapses for g^^ ^ 1. One may thus suspect that J^2u goes through zero 
simply because the approximation breaks up. In other words, the smallness of A42u in 
the QRPA could be just an artifact of the model. (One should remember that in the 
Tamm-Dancoff approximation, i.e., in the absence of the ground state correlations, A42u 
always increases with g^^.) Yet, it has been pointed out more than once that the zero of 
M 2v is not engendered by the collapse of the QRPA, but arises instead from the partial 
restoration of the S'[/(4) Wigner symmetry 0. 

It has been shown recently that within the QRPA the above behavior of the 1v am- 
plitude can be summarized as [ffl 



1 _ nVV IrpV 

M2. = M2M'' = 0) ^ \J;l , with g'o' = ^,9r^9'o', (1) 

J- ~ 9 / 9i 



where g^^ and gi^ denote respectively the zero and the pole of J^2u- Moreover, it has 
been suggested that within the QRPA the Oz/ amplitude behaves as 

Mo. = MoAJ^ = l+;g''^ = 0) ^ ' 



1 - gr>v/gl^ 

+ A^o.(J^7^1+;^7^^ = 0)(l-(?^7^7f), (2) 

where g^ ^ gi^ 0]. This means that the J'^ = 1^ component of Aiou exhibits the zero 
and the pole at the same value of g^^ as Ai2u- Thus, the theoretical estimation of Aiou, 
and therefore the determination of the limit for the effective neutrino mass < rrii, >, is 
also uncertain as that of M2u- 

Several modifications of the QRPA have been proposed in order to change the above 
behavior in a qualitative way, including higher order RPA corrections f^, nuclear defor- 
mation P, single-particle self-energy BCS terms [0 and particle number projection [0. 

^ It was found that the contributions of the odd-parity nuclear operators to the /3/32i/-decay are 
significant when compared with the experimental data [El. 
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Yet, none of these amendments inhibits the collapse of the charge-exchange QRPA. In the 
present work we show that this can be achieved by including the effect of ground state 
correlations in the QRPA equations of motion. The corresponding formalism, referred to 
as renormalized QRPA (RQRPA), was originally introduced by Rowe 0. It has been used 
recently by Catara et al., [|1^, [Tl| in the evaluation of the charge transition densities and 
properties of the charge- conserving collective states. We briefly outline below the RQRPA 
formalism for charge-exchange excitations, and discuss it within a schematic model for 

M2u- 

We begin by defining excited states \\J) that are built by the action of the charge- 
exchange operators 



^\>^J) = E [Xpn{\J)AlSJ) - V(AJ)Ap„( J) 

pn 



(3) 



on the correlated ground state |0). Here J) = and ctj^ and a]^ are quasiparticle 

creation operators for protons and neutrons. The amplitudes X and F, the eigenvalues 
ojx and |0) are obtained from the equations of motion (EM) 



(0| [5VL{\J),H,n\\J)\ |0) 
with the condition 



cja(o| \5n{xj),n\\j)\ |o). 



(4) 



^](AJ)|0) = O,for all A, J. (5) 

The usual QRPA equations result from (^) when |0) is approximated by the BCS ground 
state \BCS) and @ is ignored. In the RQRPA one takes the ground state correlations 
(GSC) introduced by in the EM (^) partially into account. First note that we have 
now 

J-i(0| [Ap„(J),4,„,(J')]° |0) = V5„„'5jj'/^p„, (6) 

with 



D 



pn 



J-i(0Mp„(J),4„(J) |0) = l-Arp-Ar„, 



(7) 



where J = \/2J + 1 and J\fp (Mn) are the proton (neutron) quasiparticle occupations 

M=jr\0\[ala,m. (8) 



The label t stands for p and n. 

We define next "renormalized" two quasiparticle operators as 



pn \ 



pn 



which satisfy the relation 



J-'{0\ Apn{J),Al,,AJ') \0)=5pp'5nn'5jJ'- 



(9) 



(10) 
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The crucial RQRPA assumption is the generahzed quasiboson approximation 



J-^ [Apn{J).^l'n'{J')\^ = J'\^\ [An(^),4'n'(^')]° |0) = ^pp' ^ J J' , (H) 

The RQRPA equations follow straightforwardly after replacing Aj,„(J) by ^],„(J) in the 
expression for V[^{J) and using (|Tl|) in the EM (|). We get in this way 



m B(j) 

B*(J) A*(J) 



X(AJ) 
Y(AJ) 



X(AJ) 
-Y(AJ) 



where 



Xp„(AJ) = Xp„(AJ)D^/2 and Yp„(AJ) = ^(AJ)^;^ 



/2 
pn ' 



(12) 
(13) 



are the renormalized amplitudes. The submatrices A(J) and B(J) are found as 



(J) = (ep + e„)5pp/5„„/ + Dp/^ [F{pn, p'n', J) {upVnUp'Vn' + VpU 



G{pn,p'n, J){UpUnUp'Un' + VpVnVp'Vn')] DI!^, , 



Bp„y„' ( J) = D^l^^ [F{pn, p'n', J) {VpUnUpiVn' + UpVnVp>Un') 



G{pn,p'n', J){UpUnVp'Vn' + VpVnUp'Un')] DJ 



1/2 
p'n' 1 



(14) 



where F and G are the usual particle-hole (PH) and particle-particle (PP) coupled two- 
particle matrix elements. 

The QRPA equations are recovered from (|T3|) and (14) by taking Dpn = 1. Within the 
RQRPA one first solves (^) in the quasiboson approximation The RQRPA ground 
state then reads 

\0) = Noe^\BCS), (15) 



with 



^5=2 E J-'[Cpnp'n'{J)Al{J)Al,AJ) 
pnp'n'J 

From (lEI) it turns out that the matrix C is the solution of 



Y^X;^i\J)Cpnp,r,iJ) = Y;„.(AJ), for all A, J. 



(16) 



(17) 



pn 



Finally, by making use of this equation one finds the quasiparticle occupations 

A/; = E ■J%'\ypn'{>^j)\'; k = y: pj-^iY,,„(Aj)p. 



XJn' 



XJp' 



The value of Dpn follows from ([^ and (P^). 
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To evaluate the transition matrix elements for the jS"^ decays 

(AJ||O(J;±)||0) = (0| [r](AJ),0(J;±)]>), (19) 

with 

0(J;±) = ^0(J;z)t±(^), (20) 

i 

we only need their two quasiparticle components 

0{J- ±) = E [Kn^J-^ ^Hn{J) + (-)"a5:( J; T)A,n{J)\ , (21) 

pn 

where 



A° (J;-) = {-yj'\r.vMO{J)\\nr. (22) 



From dl) and (^2|) one gets 

(Aj||o(j;±)||o) = jy: [a;„(J;±)x;„(aj) + (-)^a;:(J;t)y;„(aj)] z^v^. 



The corresponding total strengths are 

S{J;±) = J-'Y.\{^J\\0{J;±)m\'. (23) 

A 

Within the RQRPA the BCS equations have to be solved subject to the condition that 
|0) has on the average the correct number of particles. This requirement gives 

K = T.Jt[vt+a-2v^)Xt], (24) 
t 

Np and N„ being the number of active protons and neutrons in solving the gap equations. 

We conclude the presentation of the formalism by noting that: (a) when the factors 
Dpn, which are functions of the amplitudes Y are substituted into the renormalized ma- 
trices A and B, (|T^) becomes a nonlinear system of coupled equations for the X and Y 
amplitudes; and (b) these equations have to be solved self-consistently together with the 
new BCS conditions (0). This is the price to be paid in order to take into account the 
GSC within the QRPA problem in an appropriate way. 

We will resort now to the simplest version of the QRPA for the /5/5-decay, called the 
single mode model (SMM), in which a single RPA equation is solved with two BCS vacua 
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T^ , and only one intermediate state J'^ = 1^ enters into the play [Q. Eqs. (JT^) read in 



this case 



Apn = uJo + pppn [{ulvl + fpM^)F(pn; 1) + {ulul + vlvl)G{pn; 1) 
Bpn = 2ppPnVpUnVnUp[F{pn; 1) - G{pn; l)]£'p„, 



Dpni 



where luq = —[G{pp] 0) + G{nn; 0)]/4 is the unperturbed energy. The unbarred (barred) 
quantities indicate that the quasiparticles are defined with respect to the initial (final) 
nucleus; = Up + Vp, p~^ = u^ + v"^. All the remaining notation is self explanatory. The 
perturbed energy and Dpn are obtained by solving self-consistently the set of equations: 

u - Ia^ -52 n - 1 _ f ^yn-^ 2 _ N,/ - 3(1 - Dpn) 



pn J 



with/^3(j-2+j-2). 

The transition l3(32u matrix element is 



M 



2u 



MlDpn (^y(,^ Gipn;l)Dpn \ ^o^ ^ P^P^ W^^. , ^25) 



with M'^2u 

being the corresponding unperturbed (BCS) value. 

Numerical calculations have been performed for the ^°°Mo ^^^Ru system, where 
the appropriate intermediate state is [OgT /2{n)0gQ/2{.p)Y , and Np = 2 and N„ = 2 (Np = 4 
and N„ = 0) for the initial (final) state. We have used a 5-force (in units of MeV fm?): 
V = —A7t{vsPs + VtPt)S{r), with different strength constants Vg and Vt for the PH, PP and 
pairing channels. Thus, instead of the parameter g^^ we use here the ratio t = vf^/ tiP"*^, 
whose physical value is t = 1.5. The remaining parameters for the SMM have been taken 
to be vf = 55, vf'' = 92 and vP""'"^ = 55 [|]. The results obtained within the QRPA 
(dashed lines) and the RQRPA (solid lines) for u and for A42u are shown in Fig. 1. As 
expected, the QRPA collapses close to t = 1.5. Contrarily, in the RQRPA the energy 
decreases asymptotically when t — >■ oo. For the sake of comparison, in the same figure, 
are also presented the results for the energy of the lowest = 1"^ state and for the 
2i/ matrix element of a full QRPA calculation (dotted lines), as described in ref. 0]. 
This calculation, that involves an eleven dimensional model space, both for protons and 
neutrons, also collapses. (It is very gratifying that the simple formula ( P3| ) contains the 
main physics involved in such a relative sizable calculations.) 

In summary, we have investigated the importance of GSC effects on the solutions of 
the EM for charge-exchange excitations in the renormalized QRPA. The SMM shows that, 
contrarily to what happens in the usual QRPA, the inclusion of the GSC in the EM avoids 
collapse for physical values of the PP coupling strength. However, the amplitude A42u 



5 



still passes through zero in the RQRPA, although at somewhat higher value of t (or g^^). 
It is also evident that, in the QRPA, the physical mechanisms responsible for the zero and 
the collapse of Ai2u are not the same. The behavior of this amplitude in the RQRPA is 
not anymore delineated by Eq. (P, and the dependence of the calculated transition 
rates on g^^ is weakened. In view of Eqs. and (|^), all that was just said for the 2z/ 
mode can be extrapolated also to the Ou mode. It is well known that the contributions of 
intermediate states with J'^ 7^ 1 are quite sizeable in the neutrinoless decay for physical 
value of gP^ = 1, where it is very likely that the Aioi^lJ"^ = l"*") goes to zero even in the 
RQRPA case. But as the dynamical calculation does not collapse any more, we could now 
have more confidence in establishing the upper limit for the neutrino mass. Thus, the effect 
of the GSC in the EM appear in this context as particularly relevant. We have also found 
that a full RQRPA calculation for the A42u amplitude agrees qualitatively with the SMM 
estimate. But, in analyzing the Ikeda sum rule S{J'^ = 1+; +) — S{J'" = 1+; —) = N — Z, 
we discovered that it is not fulfilled within the RQRPA. In fact, the deviations from 
this condition grow as the GSC increase (or as the PP strength parameter increases). 
On the other hand, we have verified numerically that the similar requisite for the Fermi 
transitions is fulfilled in our formalism, when only the states J'^ = 0^ are considered in 
the Eqs. (|T8|). It should be stressed that the constraints (^4|) plays a crucial role regarding 
this point. When the usual BCS constraint on the number of particles is used |T^, the 



sum rule for the Fermi transitions is never fulfilled. That the Ikeda sum rule is necessarily 
violated in the RQRPA, when the usual BCS occupation numbers are employed, is seen 
immediately from the relation 



S{r = 1+; +) - S{r = 1+--) = W I < p| |a| |n > \\vl - v^D,^, 

pn 

which yields N — Z only when Dpn = 1. Why the Ikeda sum rule is not satisfied, even 
when the condition (^) is adopted, is still an open question. In summary, we feel that, 
before a quantitative comparison of the calculations with the experimental data could be 
done, the behavior of the sum rules in the RQRPA should be thoroughly elucidated and 
this is our next goal. 

Finally, it should be mentioned that, after our work has been completed, we have 
learned that a similar study has been performed by Toivanen and Suhonen |T3[ . 
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Figure Captions 



Figure 1: Energies uj (in MeV) of the lowest 1+ state of ^'^"Tc and the matrix elements M.2u (in 
[MeV]'"^) for the ^"^A/o "^^^Ru system. The single mode model results are indicated by the dashed 
lines for the QRPA and by the solid lines for the RQRPA. The results of a full QRPA calculation ^ are 
represented by dotted lines. 
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